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Supernova 1993 J: a veiled pulsar in a binary system? 
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The recent report of a binary companion to the progenitor of supernova 1993 J may provide 
important clues for identifying the nature of the nascent compact object. Given the estimates 
QQ ' of the progenitor mass, the potential power source is probably a pulsar rather than an accreting 

black hole. If there is a pulsar, one would expect the rotational luminosity to be stored and 
(N ■ reprocessed by the supernova remnant, but no pulsar nebula has yet been seen. The lack of 

detection of an X-ray synchrotron nebula should be taken as strong evidence against the pres- 
ence of a bright pulsar This is because absorption by the surrounding supernova gas should 
, be negligible for the light of the companion star to have been detected. A model is devel- 

' oped here for the luminosity of the pulsar nebula in SN 1993 J, which is then used to predict 

, the spin-down power of the putative pulsar If one exists, it can be providing no more than 

' ~ 6 X lO^^erg s~^. With an initial rotation period of 10-30 ms, as extrapolated from young 

galactic pulsars, the nascent neutron star can have either a weak magnetic field, B < 10^^ 
G, or one so strong, B > 10^^ G, that its spin was rapidly slowed down. The companion 
star, if bound to the neutron star, should provide ample targets for the pulsar wind to interact 



O ' and produce high-energy gamma-rays. The expected non-pulsed, GeV signal is calculated; it 

could be detected by current and future experiments provided that the pulsar wind velocity is 

^ ' similar to that of the Crab Nebula. 

■ 
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1 INTRODUCTION 

Supernova 1993J was discovered on 28 March 1993 in the spiral galaxy M81 (at a distance of 3.63 Mpc; Freedman et al. 1994). Although 
it is factor of 65 further away than SN 1987 A, it is relatively close for a SN and has been the subject of intense observational campaigns 
in a number of wavelengths regions (Chevalier 1997, and references therein). Shortly after the explosion, its spectrum was found to contain 
hydrogen lines - a type II supernova - but within a few weeks strong helium lines developed. The spectrum looked then more like a type lb 
supernova, as though the envelope of hydrogen had lifted to reveal a helium layer (Filippenko et al. 1993). This, together with the observed 
peculiarities observed in the lightcurve, led many to conclude that the dying star, identified as a non-variable red supergiant in images taken 
before the explosion (Aldering et al. 1994), must have lost a considerable amount of its outer envelope of hydrogen gas to a companion before 
it exploded (Nomoto et al. 1993; Podsiadlowski et al. 1993; Woosley et al. 1994). But no companion had been seen until recently when the 
brightness of SN 19931 dimmed sufficiently that its spectrum showed the features of a massive star superimposed on the supernova (Maund 
et al. 2004). 

Given the estimates of the progenitor mass, the remnant is probably a neutron star rather than a black hole (Nomoto et al. 1993; 
Podsiadlowski et al. 1993; Woosley et al. 1994), and the neutron star is generally expected to manifest as a pulsar. The newly born neutron 
star is characterized by its initial rotation rate and magnetic dipole moment. It will spin down, generating radiation and accelerating charged 
particles at the expense of its rotational energy. We expect the reprocessing of rotational energy to produce a bright flat-spectrum synchrotron 
nebula, but this soft emission has not been observed. 

Here we developed a simple shocked wind model for the luminosity of the pulsar nebulae in SN 1993J (Section 2), which is then 
compared with the current luminosity limits in order to place bounds on the dipole emission of the nascent neutron star and its birth properties 
(Section 4). We also investigated the impact of the soft photon field radiation of the binary companion on the central pulsar wind (Section 3), 
along with the types of observation that would help to unambiguously demonstrate whether or not the system is disrupted as a result of the 
supernova. Our conclusions are discussed in Section 4. 
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2 A MODEL FOR THE PULSAR NEBULA IN SN 1993J 

Neutron stars formed in core collapse SN explosions are thought to emit radiation as pulsars, the Crab (Rickett & Seiradakis 1982) being 
the canonical example. Pulsars are generally presumed to lose rotational energy by some combination of winds of relativistic particles and 
magnetic dipole radiation at the frequency of the neutron star's rotation. In the generic pulsar model, the field is assumed to maintain a 
steady value, and the luminosity declines as the spin rate slows down (Shapiro & Teukolsky 1983). The spindown law is then given by 
Q.{t) = no(l + t/tn)~^''^, so that Ln(t) = Ln,o(l + t/tn)~'^, where tn ^ 10^ s UiB~^^Po _3R'^'^ is the characteristic time scale for 
dipolar spindown. Bp, 15 = i3p/(10^^G) is the dipolar field strength at the poles, Re, = i?/10^cm is the radius of the light cyUnder, Po,-3 
is the initial rotation period in milliseconds, and Ln,o — 10*'' erg s'^Bp ^i^^PQ tzRt- 

The power output of the pulsar is typically assumed to be Ln, regardless of the detailed process by which the dipole emission is 
converted to the energy of charged particle acceleration and high-energy electromagnetic radiation. The pulses themselves, even when seen 
in the gamma-rays, constitute an insignificant fraction of the total energy loss. Large fraction of the power output, on the other hand, goes 
into a bubble of relativistic particles and magnetic field surroimding the pulsar (Chevalier & Fransson 1992). This bubble gains energy from 
the pulsar, and loses it in synchrotron radiation losses and in doing work on the surrounding supernova gas, sweeping it up and accelerating 
it. 



2.1 Hydrodynamic Evolution 



The density structure of the supernova gas is determined by the initial stellar structure, as modified by the explosion (Amett 1988). The 
density structure has been relatively well determined for models of SN 1987 A - it shows an inner flat profile outside of which is a steep 
power-law decrease with radius (e.g. Shigeyama & Nomoto 1990). Chevalier & Soker (1989) approximated the velocity profile by an inner 
section pi = Ar~'^t'^~^ and an outer section po = Br~"t"~^, where m = 1 and n = 9 for SN 1987 A. The two segments of the density 
profile are assumed to be continuous and they intersect at a velocity vt. For an explosion with total energy E and mass Mt, one has 
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SN 1987A was an unusual type n supernova because of its small initial radius. This is in contrast to SN 19931 where the early lightcurve 
indicated that the progenitor star was more extended, about ten times larger than that of SN 1987 A. The lack of a plateau phase in this case 
is attributed to the low mass of the hydrogen envelope, about 0.1 to 0.6 M© (Nomoto et al. 1993; Podsiadlowski et al. 1993; Woosley et al. 
1994), although the initial stellar mass was probably ~ 15M0. Woosley et al. (1994), for example, estimated that the final presupernova 
star^ had a helium and heavy element core of ~ S.TIM©, a low density hydrogen envelope of 0.2Mo, and a radius of about 4.3 x 10^^ cm. 

The core collapse of the progenitor of 1993J resulted in the deposition of about 3 x 10®° erg of kinetic energy per solar mass (Woosley 
et al. 1994). The expanded density profile shows an outer steep power-law region with n ~ 9 and an inner, relatively flat region (Blinnikov 
et al. 1998 and Fig. 4 therein). The bend in the profile occurs at a velocity of ~ 2000 km s~^. The inner density profile is clearly not a 
single power law in r. However, it is likely that hydrodynamic instabilities, similar to those believed to occur in SN 1987A (e.g. Fryxell et al. 
1991), lead to mixing and smooth out any sharp features in the density profile. Thus, a reasonable smooth, flat inner profile may be a good 
approximation. In the outer parts of the density structure, on the other hand, radiative transfer effects are more important for the explosion of 
an extended star and can lead to the formation of a dense shell in the outer layers. The shell is, however, expected out in the steep power-law 
region of the density profile, where the pulsar bubble is not likely to reach. In what follows, we assume m = 1 or 2 and n = 9. 

The first stage of evolution involves the interaction of the pulsar bubble with the inner density section. With the assumptions that the 
pulsar luminosity is constant during this phase (i.e. fage < tn) and the supernova gas is swept up into a thin shell of mass M and velocity v, 
the radius of the pulsar bubble can be written as (Chevalier & Fransson 1992) 
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km s"' for m = [1, 2], where Mt,5 = Mt/SMg, E^i = E/lQr''^ erg, Ln,39 = in/lO^^ erg s"\ and tage.ii = iage/H yr. The density 
transition velocity is vt = [5100, Q'iQQ]Ey^ M^l^'^ km for m = [1, 2], so it is likely that the shell remains within the inner part of the 
density profile. 
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The density of the uniform gas shell is p ^ [5 x IQ-'"' L^%;E;{"-M^i'/ti;:;^lA x IQ-'" Ln%; E^^''" Ml^t^^^^l] g cm"^ for 
m = [1, 2], where p oc Mf''^S"®(^"'"'/^(®"'"'L^'"''^^"'"\ The presence of the companion impUes that the supernova gas should be 



This model was derived from a 13Mq main-sequence star that lost most of its hydrogen envelope to a nearby companion (9Mq, initially 4.5 AU). 
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transparent to optical radiation (i.e. the optical depth to electron scattering should be less than unity), which gives E > W'^^ M^^^t^^^ for 
m = f and io > f 0*° E^^M'^''^ t-J^ -,^ erg s"^ for m = 2. 

For X-ray energies > 10 keV, electron scattering also provides the main opacity. However, in this case the bound electrons contribute 
as well as the free ones, so that the optical depth is always given by electron scattering irrespectively of the degree of ionization. At lower 
X-ray energies e > 3 keV, photoionization provides additional opacity (Bahcall et al. 1970) and may delay the time at which the envelope 
becomes transparent by an additional factor ~ 5(e/lkeV)~'^''^ (this estimate is for a predominantly neutral envelope of pure hydrogen). The 
delay is somewhat larger if the envelope is rich in heavy elements. The matter that is swept up by the pulsar bubble is subject to Rayleigh- 
Taylor instability and may form filaments (Vishniac 1983). The optical depth to electron scattering could then be much less if the envelope 
is sufficiently irregular that some lines of sight to the pulsar traversed relatively little envelope mass. Under such favourable circumstances, 
the opacity of the envelope is a less serious problem (except at late times when grain formation may occur in the supernova ejecta and dust 
absorption could significantly increase the opacity at visual wavelengths). 

Radio measurements suggest that any pulsar nebula in the center of SN 1993J is fainter than ~ 10"^* erg s~^ (Bietenholz et al. 2003). 
The material immediately surrounding the putative pulsar is, however, expected to be totally opaque until 



2xW^E-,'^^Mt,5T-'^%'yT formal 



at which time the radio flux (i/g = u/l GHz) will suddenly appear. Here T4 = r/10*K is the temperature of the swept-up material. The 
lack of detection of a synchrotron nebula at radio wavelengths should not be taken as strong evidence against the presence of a bright pulsar 
(Bahcall et al. 1970). On the other hand, the lack of an X-ray nebula indicates that, if one exists, its radiative power must be less than 
Lx ~2 X 10^* erg s"^ (Zimmermann & Aschenbach 2003). 



2.2 Emission Model 

A compelling model for the optical/X-ray properties of the Crab Nebula was developed by Rees & Gunn (1974). In this model, the central 
pulsar generates a highly relativistic, particle dominated wind (with Lorentz factor 7^,) that passes through a shock front and decelerates 
to match the expansion velocity set by the outer nebula. The emission from the pulsar bubble is thought to provide a larger luminosity 
source than the radiative shock front itself (Chevalier & Fransson 1992). The wind particles acquire a power-law energy spectrum of the 
form A'^(7) oc 7^^ (for 7 > 7™, where 7 is the particle Lorentz factor and 7,n is its minimum value) in the shock front and radiative 
synchrotron emission in the down stream region (Chevalier 2000). Under the basic simplification that the emitting region can be treated as a 
one zone (i.e. no spatial structure in the nebula), and assuming that a balance between injection from the shock front and synchrotron losses 
is established, the number of radiating particles at a particular 7 is given by A'^(7) — 7^"^ (7m/9B^)~^Ln7~'''"'"^' (Chevalier 2000), where 
B is the magnetic field in the emitting region, and /3 — 1.06 x 10~^^ cm'^ s~^. 

In what follows, it is assumed that the energy density in the emitting region (which is approximately determined by the shock jump 
conditions) is divided between a fraction te in particles and a fraction eg in the magnetic field. These efficiency factors are constrained 
by ee + €b = 1. With this assumption, the magnetic field in the emitting region is B = {Q^BLn/rlcY^'^ , where is the shock wave 
radius. The electron energy is radiated at its critical frequency 1^(7) — {q^B /2-Kmec), where qe and me are the electron charge and 
mass, respectively. If the electrons and positrons cool rapidly by synchrotron radiation, as thought to be the case for the Crab Nebula, 
the luminosity produced from the pulsar power is given by L„ — </>(p)ee~^eg ^^'''*7£~^rs~'''^*I/[^^'^'^*i^~''/^ (Chevalier 2000), where 
<P(P) = and 1/) = 2.8 X lO'^ in cgs units. If the particle spectrum is similar to that of the Crab Nebula (p = 2.2), 

the X-ray luminosity {u — 10^** Hz) becomes 

where cgs units are used. We note that the X-rays from a putative pulsar nebula in SN 1993J with La. ~ 10"^" erg could provide the 
additional energy input required to reproduce the observed Ha luminosity at ~350 days (Houck & Fransson 1996). 

The above relation can be used to predict the spin-down power, Lq , of the putative pulsar in wind nebula where a pulsar has not yet been 
observed (Chevalier 2000). In the case of SN 1993 J, the value Ts can be determined by the condition that the envelope should be transparent 
to optical radiation (Section 2.1). The need for a particle dominated shock suggests that ee > 0.5. We set es = te = 0.5, although there 
is little dependence to es. The value of 7^, is difficult to constraint. Here we use 7,„ ~ 3 x 10®, the value typically assumed for the Crab 
Nebula (Kennel & Coroniti 1984). A better estimate of 7™ is given in Section 3. When these parameters are substituted into equation js}, the 
XMM upper limit of Lx ^ 2 x 10^* erg s"^ (0.3-10 keV; Zimmermann & Aschenbach 2003) yields Lq < 2 x 10^" erg s"\ On the basis 
of the A5CA data, Kawai et al. (1998) found the relation log Lx = (33.42 ± 0.20) + (1.27 ± 0.17) log(Ln/10^''), where Lx is the nebular 
luminosity in the 1-10 keV range, from which we estimate Lq < 5 x lO'^'erg s~^. This estimate is in fair agreement with our predicted 
value. 
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3 THE ROLE OF BINARITY 

At the time of the explosion, the progenitor of 1993J has a mass of 5.4 Mq (with a helium-exhausted core of 5.1 A/q), the secondary has a 
mass of 22 Mq. The orbital period of the system is ~ 25 yr, and the companion has an orbital velocity of ~ 6 km s^^ (Maund et al. 2004). 
There may have been a large kick imparted by the explosion mechanism, and that would be very interesting to study, but baring that, let us 
assume the companion star is bound in an eccentric orbit with the newly born pulsar^ . The bulk of the pulsar energy Ln would be primarily 
in the form of a magnetically driven, highly-relativistic wind consisting of e~, and probably heavy ions with Lw — CLsi and C < 1- 

Under the foregoing conditions, the relativistic wind (which is likely to be undisturbed by the presence of the binary companion'^ 
since Ln 2> v'^M^) would escape the compact remnant while interacting with the soft photon field of the companion with typical energy 
df. = kT, / [nieC?) ~ 3 X 10~® (Maund et al. 2004). The scattered photons whose energy is boosted by the square of the bulk Lorentz factor 
of the magnetized wind (i.e. 6* = 27^^* ~ 6 x 1067^^ 5) propagate in a narrow 7"^ beam owing to relativistic aberration. The rate of 
energy loss of a relativistic particle moving in a radiation field with an energy density ~ L«/47ra^c is about m^c^ d'y / dt ~ —w^o-rcy^ 
in the Thompson limit (Landau & Lifshitz 1975). Here L* ~ 10^ is the luminosity of the optical star and a is the binary separation 
(a ~ 25 AU ~ 200i?, ; Maund et al. 2004). The total luminosity of scattered hard photons Lr is then equal to the total particle energy losses 
in the course of motion from the pulsar to infinity Lr = L^— ~1 — (1 + ?)^^, where 



( 3 X 106 ) ( IQSLq ) ( 200i?. ) 



denotes the efficiency in extracting energy from the relativistic outflow (Chernyakova & Illarionov 1999). The resulting radiation pressure on 
electrons in the ejecta will brake any outflow whose initial Lorentz factor exceeds some critical value 7iini < ■^{jp-Y^'^ ~ 3 x lO^Lj/gg, 
converting the excess kinetic energy into a directed beamed of scattered photons. With 7^, ~ 7iim, equation jsj yields Ln < 6 x 10"^^ erg s~^. 
As the stellar companion emits a black body spectrum, of effective temperature 8,, the local photon energy density is given by 

where zu is the soft photon energy in units of m^c^ . The scattered photons are boosted by the square of the Lorentz factor so that the local 
spectrum has a black body shape enhanced by 7^^^. As can be seen in Fig.Q the resulting spectrum is the convolution of all the locally emitted 
spectra (i.e. dn^ [r, zu\) and it is not one of a blackbody. Note that Klein-Nishina effects are important for incoming photon energies 
such that S,7„ > 1. The maximum energy of the scattered photons in this regime is 'y^rrisC^ ■ 

We note that the total luminosity emitted by the relativistic outflowing wind through the Compton-drag process in the direction of the 
observer could be highly anisotropic and may change periodically during orbital motion for an eccentric (or a highly inclined) orbit. The time 
dependence, in this case, will be very distinctive and such effects should certainly be looked for. 



4 DISCUSSION 

Despite presumptions that a neutron star may have been created when the progenitor star of SN 1993J exploded and its core collapsed, no 
pulsar has yet been seen. If one exists, its radiative power must be less than Ln ~ 6 x 10'^'"' erg or lower if the pulsar nebula has a 
high radiative efficiency. The well-known Crab pulsar and its nebula for comparison, put out 2 x 10^* erg s~^; and originally, when it was 
spinning faster, the luminosity might have been seven times greater still. Fig.|2|shows the pulsar spin-down luminosity divided by 47r times 
the square of the distance, the total pulsar energy output at Earth. The ten gamma-ray pulsars (including candidates) are shown as large 
circles. The solid curves show the dipole emission of the putative pulsar in 1993J (at t ~ 11 years) for various assumptions regarding its 
initial period and magnetic field strength. The pulsar spinning down by magnetic dipole radiation alone, with initial rotation periods of 10-30 
ms - as extrapolated for galactic young pulsars - can have a spin-down luminosity below ~ 6 x 10^^ erg (see shaded region in Fig.|2j 
only if either the magnetic field is relatively weak < 10^^ G or if it is so strong (i.e. > 10^'' G) that the pulsar luminosity decays rapidly. 
If weak magnetic moments could be ruled out in the near future by X-ray and infrared observations, then we find that, if undetected, the 
putative pulsar in SN 1993 J could be a magnetar. 

With a spectrum similar to that of the Crab (see Fig. |2j, the pulsed emission is likely to be below the detection limit of current in- 
struments and may deprive us of the opportunity to witness the emergence of a gamma-ray pulsar in the immediate future. GeV radiation 
with luminosities high enough to be detected with ARGO'^, and the Veritas^ experiment now under construction, could be produced by the 
interaction of the pulsar wind with the soft photon field of its companion provided that the system remains bound and 7^ > 10*' (i.e. similar 
or larger to than inferred for the Crab pulsar). Its detection will surely offer important clues for identifying the nature of the progenitor and 
possibly constraining whether or not kicks have played an important role. 

^ The binary is likely to remain bound after the supernova explosion provided that the kick velocity < 90 km (Brandt & Podsiadlowski 1995). 

Here and are the velocity and mass-loss rate of the stellai' companion. 
* http://argo.na.infn.it/ 
^ http://veritas.sao.arizona.edu/ 
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If there is not a neutron star in SN 1993 J, could there be a black hole instead? Had matter fallen back onto the nascent neutron star 
(with mass of 1.4 Mq) on a timescale of 100 seconds to a few hours after the explosion, enough mass may have been accreted to push the 
object over the minimum thought to be necessary for the creation of a black hole. Matter would continue to accrete onto the black hole, but 
the resulting radiation would be trapped in the outflow so that the escaping luminosity would be small. Because the cores of massive stars 
increase with initial stellar mass, this picture would be more plausible if the initial mass of the SN 1993 J progenitor star, ~15 A-Iq, were 
close to the mass limit above which stars collapse directly to black holes (e.g. Fryer 1999). 

Even if the neutron star is rotating only slowly or has a weak magnetic field, one would expect surrounding material to fall back onto 
it, giving rise to a luminosity ^ I/Edd ~ 3 x 10^® erg s~^ (see Fig. 1). The remaining possibility of a weak pulsar with little surrounding 
mass will be difficult to rule out. The thermal emission from a newly formed neutron star gives a luminosity of ~ 3 x 10^* erg s^^. The 
neutron-star emission should have a characteristic soft X-ray emission, but its detection will be arduous. If future observations fail to identify 
a pulsar in 1993J, then the youngest pulsar that we know will still be PSR J0205-F6449 (Camilo et al. 2002), associated with SN 1181. 
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Figure 1. Representative spectra produced by the interaction of tlie pulsar wind with the radiation of the stellar companion. The luminosity of the gamma-ray 
emission Lp is shown for 7^, = 10^, 10"^, 10'', lO'' (dashed Unes). The spectra of the Crab and Vela pulsars are shown for comparison (adapted from 
Thomson 2003). 
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Figure 2. Pulsai' observability as measured by the spin-down energy seen at Earth. Small dots represent pulsars with no gamma-ray emission, while large 
filled symbols show all known gamma-ray pulsars (including candidates; Thompson 2003). The power-output of the pulsar wind, at the distance and age (here 
assumed to be 11 years) of SN 1993J, is shown for B = lO", lO", lo"^ lo" and iQib q ^f^^^ ^g^t) 

as a function of the assumed initial 

period. The putative pulsar in SN 1993 J, spinning down by magnetic dipole radiation alone, can have a wind nebula X-ray luirtinosity below the XMM limit if 

Lq < 6 X 10^^ erg (shaded region). 



